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ON : ABSTRACT 

Q^, Calculations of the rate of local Primordial Black Hole explosions often assume 



that the PBH's can be highly concentrated into galaxies, thereby weakening the 
Page-Hawking limit on the cosmological density of PBH's. But if the PBH's are 
concentrated by a factor exceeding c/{HoRo) « 4 x 10^, where i?o = 8.5 kpc is 
the scale of the Milky Way, then the steady emission from the PBH's in the halo 
^ ■ will produce an anisotropic high latitude diffuse gamma ray intensity larger than 

. the observed anisotropy. This provides a limit on the rate-density of evaporating 

, PBH's of ^ 0.4 pc~^yr^^ which is more than 6 orders of magnitude lower than recent 

[ experimental limits. However, the weak observed anisotropic high latitude diffuse 

l/^ \ gamma ray intensity is consistent with the idea that the dark matter that closes the 

' Universe is Planck mass remnants of evaporated black holes. 

O ' Subject headings: gamma rays: observations - cosmology: dark matter - Galaxy: halo 

1. Introduction 

X 

\—i ' The average density of primordial black hole's (PBH's) in the Universe is constrained by 



the Page-Hawking (1976) limit on the diffuse gamma ray intensity, since the Hawking (1974) 
radiation from PBH's produces copious gamma rays. Halzen et al. (1991) computed the photon 
spectrum from a uniform density of PBH's with the initial mass function ni{mi)dmi oc m~'^'^dmi, 
and found that QpBuh'^ < (7.6 ± 2.6) x 10~^ with h = i?o/100 km/sec/Mpc, or an average mass 
density of PBH's of ppbh < 1-43 x 10^"^^ gm cm^'^, and an average number density of PBH's of 
N = 10"^ pc~3_ 

Now I wish to calculate the maximum allowed density of PBH's in the halo of the Milky Way. 
The halo mass density is given by 

Ph = = 8.4 X 10-25 gm cm-3 (1) 

at the position of the solar circle, Ro = 8.5 kpc, for a circular velocity of Vc = 220 km/sec. Since 
the halo density is (4.5 x l{f')/Q.h? times higher than the average density of the Universe, one 
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could hope that the PBH's would have a higher density in the halo of the Milky Way which 
would make it easier to detect the explosions caused by their final evaporation in high energy 
gamma ray experiments. In fact, Halzen et al. (1991) considered concentration factors up to 
C = Ph/p = l-36//i^ X 10^ by assuming that the PBH's were as highly concentrated as the luminous 
matter. Cline & Hong (1992) considered local densities as high as N = 10^'^ pc^^ by assuming 
that some of the gamma-ray bursts observed by BATSE were PBH explosions. 



2. Calculation of Anisotropy 

However, the PBH's in the halo will contribute an anisotropic diffuse gamma ray intensity 
that will be much easier to measure for instrumental reasons than the isotropic intensity. A rough 
order of magnitude estimate for this anisotropic signal is a fraction (HqRo/c times the isotropic 
background. In order to improve this calculation, I need to compute the local average cmissivity 
from the Page-Hawking limit which is integrated over time. For simplicity, I will do this calculation 
for the bolometric gamma ray intensity, though a frequency-dependent calculation would give a 
better limit. The mass spectrum of PBH's produced by a scale invariant perturbation spectrum 
in the early Universe is ni{mi)dmi oc m^'^'^dmi, and the initial mass of a PBH that is evaporating 
at time i is mj oc 

^V3. The total comoving density p/(l -|- z)^ of PBH's scales like 

p/il + zf oc j^^^^mm~^-^dm oc t~^/^ (2) 

and thus the comoving luminosity density scales like t"^/^ oc (1 + zY^^ for U = 1. The bolometric 
intensity is related to the comoving emissivity by 

, f jCMjt) f Jo(l+^)^/^ C 

This gives the relationship between the current average emissivity and the isotropic integrated 
intensity, 

^iso — ■ (4) 

Now I will calculate the emission from the halo of the Milky Way. If the emissivity of the halo is 
jh, then its contribution to the anisotropic intensity at angle 6 with respect to the Galactic center 
is 

laniso = / 3h{R)ds = -^—^jh{Ro)Ro (5) 
J Sill u 

for a spherical halo with density following the singular isothermal sphere model: p oc r~^. This is 
a special case of an "isothermal" halo with core radius rc, flattening q, and density varying like: 

p= M + ^c) (6) 



EGRET E > 100 MeV photons/m^/sec/sr 
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Fig. 1.— Gamma rays with E > 100 MeV from EGRET. 

with R and z being cylindrical coordinates. This gives an anisotropic intensity of 



la 

with 



j jhds = jh{Ro)Ror){l,b,rc/Ro,q) 



nil, rJR.,,) = (l + C-c/^)^)f;^/2 + 'a^-(c°B>'/sin>")] 

sin^'Yl + - \)sm^h 



where cos 6' = (cos / cos b)/ \J\ + [q ^ — 1) sin^ h and sin 6" = v^l — cos"^ 0' + r'^/R^. 

Most of the uncertainty in the endpoint of primordial black hole evaporation cancels out 
the ratio of laniso/hso which is given by 

-^aniso /i i_ i rt \jh\-^o) '^HqRq 

— =r](l,h, Tc/Ro , q) : r 

^iso Jo 4C 

3. Comparison with Data 



Using the combined Phase I EGRET sky maps in the directory 
/compton_data/egret/combined_data on legacy.gsfc.nasa.gov, I have constructed the rate 
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map in Figure |T|. This shows the rate of gamma rays with E > 100 MeV. Figure || has been 
smoothed with a 2° FWHM quasi-Gaussian test function (Wright et al. 1994), but the fits are 
based on unsmoothed maps. 

PBH's are not the only, or even the primary, source of galactic gamma ray emission. The 
process 

PCR+PISM ^ V+V + '^° P + V + '^l (10) 

where a cosmic ray proton pcR hits an interstellar medium proton pisM is the dominant source 
of galactic gamma rays. Digel, Hunter & Mukherjee (1995) find an average emissivity of 
(1.7 =b 0.1) X 10~^^ photons/sec/sr /proton for E-^ > 100 MeV in the Orion region which is similar 
to interstellar medium emission rates elsewhere in the solar neighborhood. While the cosmic ray 
density appears to be higher in the inner galaxy, leading to a higher emissivity per ISM proton, I 
can avoid this complication by not using regions close to the galactic plane in my fits. I have used 
the 100 /xm intensity as a proxy for the column density of ISM protons. This will automatically 
include the II2 and H II components of the ISM that are not measured by the 21 cm neutral 
hydrogen line. The 100 //m emissivity per proton will depend on the local dust properties and 
radiation field, but avoidance of the galactic plane minimizes these complications as well. The 100 
^m intensity at high \h\ is approximately 0.6 x 10^^'^ MJy/sr/proton/cm^, so one expects about 
0.03 photons/m^/sr with E > 100 MeV for each MJy/sr of 100 /im intensity. 

When I fit the gamma ray map to the form 

/-,(/, h) = jh{Ro)Ro r]{l, b, rjRo,q) + ^hooH, h) + esc |6| + h + e{l, b) (11) 

diioo acsc|o| 

with /loo being the 100 fiui intensity from the DIRBE instrument on COBE with a model of 
the zodiacal light removed, and e being the residual, I get the results shown in Table |l|. Missing 
parameter values are fixed at their default values, which are zero except for Vc/Ro = 0.01 and 
g = 1 (which approximate the simple singular isothermal sphere model.) The isotropic intensity 
is required to be non-negative. These fits were done using the least sum of absolute errors 
X] |e| instead of least squares fitting to avoid the effect of bright sources, and an elliptical region 
surrounding the Galactic Center with h = ±30° and / = ±45° was excluded along with \b\ < 14.5°. 
The excluded region is marked on Figure |^. The pixels used for this fit were COBE DMR 
pixels with 4328 pixels outside the exclusion region, so differences in ^ |e| greater than 0.1% are 
statistically significant. Simulating fits to 100 skies based on the best fit model in Table |l| along 
with the observed residuals gave values of jh{Ro)Ro with a standard deviation of 0.00165, which 
is negligible when compared to the range ol jh{Ro)Ro obtained using different galactic tracers and 
halo models. 

The coefficient found for the /loo term is consistent with the result expected from studies 
of the interstellar medium in the solar neighborhood (Digel et al. 1995), but when csc|6| or 
flattened haloes are allowed to absorb some of the galactic flux the coefficient is slightly lower but 
not unreasonable. Fits that do not include esc 16! as a tracer tend to favor flattened haloes with 
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large core radii, which makes the halo model ij approximately proportional to esc Models that 
do include a separate esc \b\ term favor spherical, singular isothermal sphere haloes. Because a 
flattened halo has a smaller thickness at high the fitted halo density is higher for the flattened 
models. The large core radius in the flattened models causes the halo to produce a large monopole 
contribution to the intensity, so the isotropic term goes to zero in the flattened models. 

Models similar to the large core radius spherical haloes proposed for gamma ray burst (GRB) 
sources do not provide an anisotropic intensity, for the simple reason that the GRB's are observed 
to be isotropic. In these models the halo density is limited by the isotropic intensity, and scales 
like ]h oc i"c^- Thus large core radius models are not signiflcant for placing an upper limit on the 
local halo density. 

Different methods of fltting for the galaxy give very different estimates for the isotropic 
background, but the estimate for the local halo density jh{Ro) is slightly more stable. Thus it is 
more appropriate to normalize to the Halzen et al. )1991) model upon which their calculation of 
the Page-Hawking limit is based instead of the uncertain isotropic background. Integrating Figure 
3 of Halzen et al. (1991) for E > 100 MeV gives a flux of 0.06 photons/m^/sec/sr which is in the 
range of the isotropic intensities from the fits in Table |l|, and is also consistent with the darkest 
sky intensity of 0.15 photons/m^/sec/sr given by Bertsch et al. (1995). 



Comparing these results with Equation ^ implies that 

jh{Ro)Ro ^ 0.023-0.15 
(4/3)jo(c/i/o) ~ 0.06 

or 



0.4-2.5 (12) 



C = X 0.4-2.5 = (2-12)//i X 10^ (13) 

Using this C I get a local density of PBH's of iV,, = C x 10^ = (2-12)//i x 10^ pc^^. 

An alternative fit is shown in Figure ^. The data with |6| > 19.3° was binned into twenty bins 
equally spaced in cos^. Within each bin, the mid-average, or average of the two middle quartiles 
of the data sorted by intensity, was taken. The filled points in Figure ^ are these mid-averaged 
intensities. Since large | cos^| only occurs for small |6|, the mid- average value of esc \b\ was also 
computed in each 9 bin. The mid-averaged intensities were then flt to the form 

I^ = Io + CSC |6| + jh{Ro)Ror]{l, b, 0, 0) (14) 

giving coefficients lo = —0.005 photons/m^/sec/sr, a slope of dl-y/dcsc\b\ = 

0.14 photons/m^/sec/sr, and a local halo emissivity oi jh{Ro)Ro = 0.037 photons/m^/sec/sr. This 

fit is the curve in Figure and the contribution of the esc |6| term is shown as the open circles. 
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Table 1. Halo fits to the EGRET map. 



jh{Ro)Ro 


dIj/dliQo 


dij/dcsc \b\ 


rc/Ro 




lo 


E|e| 












0.2818 


413.6 


0.0454 










0.2006 


393.5 


0.1504 






2.03 


7.5 





303.2 






0.1075 






0.1005 


319.7 


0.0409 




0.1024 






0.0364 


298.9 


0.0512 




0.0952 


0.01 


1.5 


0.0413 


298.8 




0.0330 








0.1692 


310.7 


0.0229 


0.0317 








0.1345 


303.8 


0.0743 


0.0183 




2.89 


5.0 





280.9 




0.0218 


0.0574 






0.1099 


289.4 


0.0306 


0.0173 


0.0656 






0.0578 


278.1 



Note. — I-y, /o, j/i-Ro and e are in photons/m^/sec/sr, while /loo is in 
MJy/sr. 
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Fig. 2. — Fit of a esc |6| plus halo model to the E > 100 MeV intensity. 9 is the angle between the 
line-of-sight and the galactic center. 



4. Discussion 

If PBH's were strongly concentrated in the halo of the Milky Way, they would produce a 
large anisotropic gamma ray flux which could easily be observed. A weak anisotropic signal of the 
predicted form is present. With the resulting value for the concentration factor C, ~ (2-12)//i x 10^, 
I can use the Halzen et al. (1991) calculation of the average PBH explosion rate to estimate that 
the explosion rate of halo PBH's is 



dn 
lit 



< 0.07-0.42 pc'V"^- 



(15) 



Since there are several possible emission mechanisms other than PBH's which could be located in 
a galactic halo, I have taken the rate from the fits as an upper limit. The dominant uncertainty 
in this limit is the physical thickness of the galactic density enhancement, and this is reflected in 
the range of models considered in Table |l|. The highest densities correspond to flattened haloes, 
and for a collisionless species like PBH's the highest likely flattening corresponds to an E7 galaxy 
shape with = 11. The best fit values obtained here when not including a separate esc |6| term. 
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q ^ = 5 or 7.5, are both equivalent to an E6 galaxy shape. The observed flattening of the dark 
matter in polar ring galaxies is in this range (Sackett et al. 1994). 

The uncertainties in modeling the last few seconds of PBH lifetime have very little effect 
on the limit derived here, because the diffuse gamma ray flux comes primarily from PBH's with 
masses M > 10^^ gm and temperatures kT < 100 McV, and radiation under these conditions is 
well understood. However, the behavior of PBH's at higher temperatures is not so well known, 
and different models of the final burst can give vastly different detection limits. For example, 
the recent EGRET limit of dn/dt < 0.05 pc~^yr~^ (Fichtel et al. 1994) assumed that the last 
6 X 10^^ grams of PBH rest mass evaporate producing 10^^ ergs of 100 MeV gamma rays in less 
than a microsecond based on the Hagedorn (1970) model for high energy particles, while in the 
standard model of high energy particle physics it takes > 10^ years to evaporate this mass (Halzen 
et al. 1991). The particular technique used in Fichtel et al. (1994) would not be sensitive to 
standard model bursts. Porter &; Weekes (1978) derived a limit of dn/dt < 0.04 pc~^yr~^ using 
the Hagedorn model, but this limit is weakened to dn/dt < 5 x 10^ pc~^yr~^ in the standard 
model (Alexandreas et al. 1993). Cline & Hong (1992) proposed that an expanding fireball could 
convert much of the 10^^ ergs from a Hagedorn model burst into MeV gamma rays detectable by 
the BATSE experiment. Given the limit on dn/dt above, BATSE would have to be able to detect 
PBH explosions out to distances 0(10^^) cm if even a few percent of the BATSE bursts were due 
to PBH's. Other sensitive limits on the PBH explosion rate density (Phinney & Taylor 1979) 
depend on conversion of the last 10^^ grams of PBH rest mass into an expanding fireball that 
produces a GHz radio pulse by displacing the interstellar magnetic field (Rees 1977). However, in 
the standard model it takes a few days for the last 10^^ grams to evaporate (Halzen et al. 1991), 
so no radio pulse is generated. 

Similarly, variations in the initial mass function of PBH's do not affect the ratio of the local 
emissivity to the rate density of evaporating bursts, because both the 10^^ gm PBH's radiating 
the diffuse gamma rays and the 10^ gm evaporating PBH's were initially formed with very similar 
masses near 10^^ gm. Thus changing the slope of the initial mass function has very little effect 
on the ratio of their abundances, even though such a slope change has a large effect on the 
Page-Hawking limit on Q,pBH- 

The limit derived in the paper on the local rate-density of evaporating PBH's provides a very 
difficult target for all techniques to directly detect standard model PBH explosions. For example, 
the CYGNUS experiment presented a limit dn/dt < 8.5 x 10^ pc~^yr~^ (Alexandreas et al. 1993), 
and my estimate is 6-7 orders of magnitude smaller. 

On the other hand, the ratio of the halo density of PBH's to the Page-Hawking limit is quite 
close the ratio of the total halo density to the critical density in the Universe. This suggests that 
PBH's could be tracers of the Cold Dark Matter (CDM). This would naturally occur if the PBH's 
were the CDM, but this requires either a modified mass spectrum with an enhanced abundance 
of PBH's with M > 10^^ gm, or else that evaporating PBH's leave behind a stable Planck mass 
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remnant (MacGibbon 1987; Carr, Gilbert & Lidsey 1994). In general, the fits in this paper give 
concentration ratios that are slightly too high for this hypothesis if = 1 and h = 0.8. The lowest 
ratio of jh{Ro)Ro to Iq in Table [l| gives a concentration of Clow = S/h x 10^, while if PBH's trace 
the dark matter the expected ratio is (4.5 x 10^)/O/i^. These are equal only if il/i = O.BAQow/C- 
For a more typical low value of C ~ 2//i x 10^ (corresponding to dn/dt « 0.1 pc^^yr^^), the 
required value of fi/i is ~ 0.2, which is consistent with theories of large-scale structure formation in 
CDM (Peacock & Dodds 1994) or CDM+A (Efstathiou, Sutherland k Maddox 1990) models. If 
this admittedly weak correspondence is correct, then 100 MeV gamma rays are providing the first 
non-gravitational evidence for CDM. Any model for dark matter that gives a 100 MeV gamma ray 
emissivity proportional to the density, such as particles with a very slow radiative decay, would 
also be supported by this correspondence, while models with emissivity proportional to p^, such 
as annihilating particles, would not. More data and better galactic modeling are needed to test 
this exciting possibility. 

One possible test is to look for gamma rays from large concentrations of dark matter. The 
flux from the Galaxy, F = J I cos Odd, in Figure Q is F = 1.1 7/m^/sec while the flux from the 
halo term of the best fitting model in Table [l| is -F = 0.38 7/m^/sec. Scaling the latter flux to the 
mass and distance of M87 (Stewart et al. 1984), assuming a constant gamma ray luminosity to 
mass ratio, predicts a flux of 4 x 10~^ 7/m^/sec for E > 100 MeV which is only 10 times lower 
than the limit of 40 x 10~^ 7/m^/sec reported by Sreekumar et al. (1994). Detection of a gamma 
ray flux from clusters of galaxies that correlates with dark matter column density instead of gas 
column density would support of association of dark matter and PBH's. 

This research has made use of data obtained through the Compton Observatory Science 
Support Center GOF account, provided by the NASA-Goddard Space Flight Center. 
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